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We investigate the chiral magnetic effect (CME) under a strong magnetic field B = BqXs at low 
temperature T < T*. For this purpose, we employ the instanton vacuum configuration with the 
finite instanton-number fluctuation A, which relates to the nontrivial topological charge Q t . We 
compute the vacuum expectation values of the local chiral density (p x ), chiral charge density (n x ) 
and induced electromagnetic current (j^), which signal the CME, as functions of T and Bo- We 
observed that the longitudinal EM current is much larger than the transverse one, |j±/j||| oc Qt, 
and the (n x ) equals to the | (j'3,4) | - It also turns out that the CME becomes insensitive to the 
magnetic field as T increases, according to the decreasing instanton, i.e. tunneling effect. Within 
our framework, the instanton contribution to the CME becomes almost negligible beyond T ss 300 
MeV. 

PACS numbers: 12.38.Lg, 14.40.Aq 

Keywords: chiral magnetic effect, topological charge, P- and CP-violations, instanton vacuum configuration 



I. INTRODUCTION 



Investigations on the nontrivial QCD vacuum structure have been one of the most important and dedicated subjects 
in modern particle and hadron physics for several decades. We note that the QCD ground state consists of infinitely 
degenerated vacua (0- vacuum), and each degenerated vacuum can be enumerated by its own integer topological 
number, i.e. the Chern-Simon number ncs- According to the Atiyah-Singer theorem, two different vacua can be 
related to each other by a nontrivial topological charge Q t = ncs(^ = 00 ) — n cs(t = — °°); which stands for 6 =/= 
and signals the breakdown of the P- and CP-invariance of the vacuum. Since each vacuum is separated by a certain 
potential barrier, a nontrivial Q t means quantum mechanical tunneling or passing over the barrier with enough 
energy. These novel mechanisms are known as instantons [H |2] and sphalerons [3j H] , respectively. As temperature 
(T) increases, the sphaleron contribution prevails over that of the instanton, since the tunneling rate is diminished 
as oc exp(— T 2 ) by the screening [2J. However, since the instanton contribution remains still considerable below the 
chiral restoration T, w 200 MeV, which can be lowered more with the inclusion of dynamical quarks [5j|6j[7], one 
may expect a sizable instanton contribution for the nontrivial Qt at low T. 

Now, we want to mention the recent discussions on how to measure this interesting QCD vacuum property, demon- 
strated by the nontrivial Q t , in experiments. In Refs. |8J, it was suggested that the nontrivial Qt can be probed by 
measuring asymmetric electric currents, due to an event-by-event P- and CP-violation, in the non-central heavy-ion 
collision experiments such as RHIC, FAIR, and LHC. Their main idea is that unequal numbers of the chirally left- 
and right-handed quarks, Nl — ATr. 7^ proportional to the Qt in terms of the axial Ward-Takahashi identity [5], 
can produce the asymmetric electric current under a strong magnetic field B, which is generated by the collision 
and perpendicular to the reaction plane. Interestingly enough, a possible experimental evidence, which is the electric 
charge separation, was reported recently by STAR collaboration at RHIC [TO] . Moreover, in Ref. [TT], using the 
tadpole-improved SU(2) quenched lattice simulation, the electromagnetic (EM) current fluctuation as well as local 
chirality were explored at finite T, resulting in that the fluctuations become insensitive to the external magnetic field 
as T increases. Moreover, they supported an instanton-like gauge configuration for the CME at low T, since the 
longitudinal component of the j M is drastically enhanced by it. In Ref. 12J, it was also shown in four different ways 
that the electric current, induced by the external magnetic field with the nontrivial Qt, reads 

j = ~^B, (1) 

where the fi x represents the chiral chemical potential and contains the information on the Q t inside of it. Note that 
we set the quark electric charge to be unity in Eq. ([IJ for simplicity. In fact, this is a consequence of the EM axial 
anomaly [IB] , 
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In the present work, as mentioned above, we investigate the CME at low T, employing the instanton vacuum 
configuration with the nontrivial Q t , which relates to the finite instanton-number fluctuation, i.e. the number dif- 
ference between the instantons and anti-instantons in the grand canonical ensemble, A = Nj — Nj ^ resulting in 
the CP-violation. To this end, we first write the effective action derived from the instanton vacuum configuration 
as a functional of A [14] . Using generic functional external-source and linear Schwinger methods, we compute the 
relevant physical quantities, the vacuum expectation values of the local chiral density (p x ), chiral charge density 
(n x ), EM current (ju), induced by the external magnetic field, as the signals of the CME. In order to consider the 
T-dependence of the relevant quantities, we also employ the Harrington-Shepard caloron to obtain the T-dependent 
constituent-quark mass. Fianlly, the standard fermionic Matsubara formula is taken into account for the anti-perioidc 
sum over the Euclidean time, i.e. T. We consider all the relevant physical quantities up to C(A 2 ) ~ 0(Q 2 ). 

As a result, we can obtain the correct expression for the induced EM current as given in Eq. (JT|) and observe 
that the longitudinal component of the EM current is much larger than the transverse one, satisfying the relation 
Ij'-l/jhI oc Qt- Moreover, it turns out that the (n x ) equals to the |(j3,4)| as long as the Lorentz invariance remains 
unharmed. These results are compatible and consistent with those given in Refs. [HI [12]. We also find that the CME 
becomes insensitive to the magnetic field as T increases, according to the decreasing instanton, i.e. tunneling effect, 
as observed in the lattice QCD simulation [11) . Within our low-T instanton framework, the instanton contribution to 
the CME becomes almost negligible beyond T w 300 MeV. 

We organize the present work as follows: In Section II, we make a brief discussion on the effective action derived 
from the instanton vacuum configuration with the finite instanton-number fluctuation A ^ 0. The relevant physical 
quantities, which signal the CME, are defined and computed with the effective action in Section III. All the ingredients 
obtained in the last Sections are extended to the case at finite T in Section IV using the Matsubara formula and the 
Harrington-Shepard caloron. In Section V, numerical results and related discussions are given. Finally, Section VI is 
devoted to summary and conclusion. 



II. EFFECTIVE ACTION FROM THE INSTANTON VACUUM 



In this Section, we briefly introduce a P- and CP- violating effective action 5 ft, derived by Diakonov et al. from 
the instanton vacuum configuration in the large iV c limit at zero temperature (T = 0) [14] . Employing a dilute grand 



canonical ensemble of the (anti)instantons with a finite instanton-number fluctuation, A 
corresponds to a CP-violating vacuum, but a fixed total number of the pseudo-particles 
be written in momentum space with Euclidean metric as follows: 
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where we have used Nf = 1. However, the extension to an arbitrary Nf is just straightforward. The C stands for an 
irrelevant constant for further investigations, whereas the N±/V for the (anti)instanton packing fraction proportional 
to the inverse of the average inter-(anti)instanton distance 1/P 4 « (200 MeV) 4 . The A± denotes a Lagrangian 
multiplier, which was employed to exponentiate the 2iV/-'t Hooft interaction in the effective action pp. The average 
instanton size in the dilute instanton ensemble is assigned as 1/p ~ 600 MeV, while the m indicates a small but finite 
current-quark mass for the SU(2) light-flavor sector (m — > 0). The F(k) denotes the quark form factor originating 
from the non-local quark-instanton interactions and is defined as 



F(k) = 
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(3) 



where the /„ and K n stand for the modified Bessel functions. We, however, will employ a parameterization of this 
from factor for convenience in the numerical calculations. 

From the effective action, we can obtain the following two self-consistent (saddle-point) equations with respect to 
the A±: 
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where the A± is approximated as Mo(l ± S) in the last line of Eq. |4]) with account of the fact A <C N in the 
thermodynamic limit [14] . The momentum-dependent constituent-quark mass is defined as M(k) = MQF 2 (k) pQ. By 
adding and subtracting the instanton (+) and anti-instanton (— ) contributions in Eq. Q, we arrive at 

N mM N c f d 4 k (1 + S 2 )M 2 

V 2ir 2 p 2 c J ^Tr^tf + il + S^M 2 ' K > 
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Taking into account that 5 -C 1, A -C N, and using Eq. ([6]), we can obtain an expression for S as a function of relevant 
parameters: 



/ 2n 2 p 2 \ A 
\mM N c ) V 



(7) 



This equation tells us that the 6 contains the information on the instanton-number fluctuation A at a certain scale p, 
which is about 600 MeV in the present framework. Taking into account all the ingredients discussed so far, finally, 
we can write the relevant effective action with A ^ for further investigations: 
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where the Tr Ci / )7 denotes the trace over color, flavor and Lorentz indices. 

Now, we are in a position to discuss the relation between the nontrivial topological charge Qt, as a source of 
the CME, and the instanton number fluctuation A. According to the axial Ward-Takahashi identity [9], the Qt is 
proportional to the number difference between the chirally left- and right-handed quarks, Q t <x Nr — A^^. Hence, 
the nontrivial Qt indicates the chirality flip. Note that, similarly, if a chirally left-handed quark is scattered from 
an instanton to an anti-instanton, the quark helicity is flipped to the right-handed one, and vice versa. This means 
that the nonzero A results in Nr — Nl ^ 0. In this way, the Qt can be considered to be proportional to the A: 
Qt ~ A [TJ[T5]. As a consequence, we can study the CME using the effective action in Eq. Q as a functional of Q t , 
more explicitly S cx A. The S in Eq. Q can be rewritten in terms of a real and small parameter e, which satisfies the 
condition lei < 1, for convenience as follows: 



\mM N c 



eN 
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Using the typical values for the parameters for the flavor SU(2) instanton model listed in Table [I] one can obtain 
5 « 13.57e. 



R 


P 


(N/V) 1/4 


M 


m 


1 fm 


1/3 fm 


200 MeV 


350 MeV 


5 MeV 



TABLE I: Instanton parameters, constituent and current quark masses in vacuum. 



III. THE CHIRAL MAGNETIC EFFECT IN THE MAGNETIC FIELD 



In this Section, we compute relevant physical quantities for the CME in the presence of the external magnetic field, 
using the effective action given in the last Section. They are the vacuum expectation values (v.e.v.) of the local chiral 
density p x , chiral charge density n x , and electromagnetic (EM) current j^, which are defined as follows: 

p x (x) = iq^(x)j 5 q(x), n x (x) = iq^(x)j 4 j 5 q(x), j^(x) = iq* (x^^x). (10) 

The p x and n x represent the strength of the parity breaking in a system, whereas the j M here the EM current induced 
by the external magnetic field. It is very convenient to compute these quantities employing the effective action with 
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corresponding pseudoscalar (V), axial vector (A), and vector (V) external sources. Hence, the effective action can be 
rewritten as: 
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We will take N c = 3 and Nf = 2 throughout the present work. The indicates the covariant quark momentum, 
gauged by the photon field as k^ + A^, in which we set the quark electric charge unity for convenience. Note that we 
have written the external sources V and A±, which correspond to the p x and n x , to be proportional to the 5, since 
they signal the nontrivial topological charge Q t as the parity-breaking quark mass Sj^M does. 

First, we calculate the v.e.v. of the chiral density, using the standard functional method as follows: 
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where we have performed the trace of over the color and flavor indices. Since we are interested in the response of 
the nonperturbative vacuum to the external magnetic field for the nonzero topological charge Qt, we only collect the 
terms proportional to the 5 and the photon field strength tensor F^. In order to perform the trace over the Lorentz 
index in the r.h.s. of Eq. ( 12 1 under the external EM field, we employ the linear Schwinger method [HI [19]. According 



to that, the quark propagator in the presence of the instanton background can be written up to 0(F^ U ), which is 
equivalent up to O(Bq), as follows: 
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Here, a ■ F = UavFuv- We define relevant functions related to the momentum-dependent quark mass: 
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Note that these functions are appropriate parameterizations from Eq. ^ as mentioned in the last Section. After a 
straightforward trace manipulation, one arrives at a simple expression for the chiral density: 

(p x ) F ,s = 25 2 T a Bl (15) 

where the subscripts F and 5 in the l.h.s. of Eq. |i"5"| ) indicate the fact that we picked up only the terms proportional 
to the Fpv and 5 as explained above. In deriving Eq. ( |15| , to make the problem easy, we assumed a static external 
EM field, which is assigned for example as in Ref. [TT] : 
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where the external EM field consists of the classical (A ^) in the symmetry gauge and fluctuation (A^ uc ) parts. We 
take the 01^,4 as small and constant EM potentials, satisfying <9 M ai^4 = 0. In this assignment for the A^ 7 we have 
only the constant magnetic field in the spatial X3-direction: B = B Q x^- Thus, among the EM field strength tensors, 
F12 = Bo and its dual remain finite, while others disappear. Note t hat we define a dual field strength tensor in 
Euclidean space as F^ v — ^e^upo-Fpa- The relevant integral T a in Eq. (151 reads: 



T a = —AN c Nf 



d A k M (MM + |) 
(2tt) 4 [k 2 + (1 + 5 2 )M 2 } 2 



(17) 



As understood from Eq. ( 15 1, the parity breaking of the vacuum becomes enhanced quadratically with respect to the 
S and Bq. 

Similarly to the p x , we can compute the v.e.v. of the chiral charge density using the following functional derivative 
with respect to the external axial-vector field: 
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By solving Eq. (18) and picking up relevant terms, one is lead to a compact expression for it as follows: 

(ti x )f,s = 5FbtAv P a{iA<j)F vp = 2SJ 7 b (iA 3 )B , 



(19) 



where we have used the fact that only the F34 = Bo is nonzero in our symmetric-gauge EM field. From Eq. (19), we 
observe that the chiral charge density increases linearly with respect to the 8 as well as the Bo. The integral T b is 
written as: 
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(20) 



Finally, we attempt to compute the v.e.v. of the EM current, induced by the external magnetic field in the presence 
of the nontrivial topological charge: 
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Expectedly, we can obtain a very similar expression for the (ju)F 5 to that for the v.e.v. of the chiral charge density: 



(7m>f,« = -[3^^ + 25F^]{iA v )T b , 



(22) 



where the definition of the T b is the same with that given in Eq. (|20j) . Considering Eq. ( 22 1 and assuming finite 
values for A\^,±, we can write the induced EM currents for Q t ^ separately for its transverse (_L) and longitudinal 
(||) components: 



0'i,2>F,a = j± = T3(5 2 ^ b (a 2il )B , 
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(23) 



The jj| increases linearly as the B grows, which is the indication of the CME as in Eq. (jTJ). Equating the transverse 
and longitudinal components, we have the ratio as follows: 
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0), in other words, the strengths of 



From Eq. (24), if we assume that | j4_i,2 | ~ 1^-3, 4I along the i^-direction (x\ 
the fluctuations 01^4 in Eq. ( 16 1 are almost the same to each other, it is obvious that the transverse components of 
the induced EM current are much smaller than those of the longitudinal ones, \jx/j\\ I *C 1, due to the fact that 5 <C 1. 
This tendency is just consistent with that given in the recent lattice QCD simulation [11 , in which the instanton-likc 
contribution turned out to be crucial for this large difference between the longitudinal and transverse components of 
the EM current for the nontrivial Qt. 

Considering the general expression for the induced EM current in Eq. (JT|), the induced EM current along the 
X3-direction in Eq. ( 23 1 can be rewritten as 
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1/3 fm and M « 0.35 GeV to the T b for S 
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(26) 



in Eq. ( 20 1 , Eq. ( 26 1 can be estimated as 



Substituting p 

fi x w 0.92 x (iSAi). As can be seen from the effective action in Eq. (Ill, the term 8A4 plays the role of a chemical 
potential. Hence, this estimation shows a correct behavior of the 5A~[: \i x — iSA^. About 8% shortage may be 
understood by the fact that there can be additional contributions to Eq. (26 1, if we expand the propagator in Eq. (13 1 
more to higher orders in 5. Moreover, comparing Eqs. (19 1 and (23 1, it can be easily shown that the v.e.v. of the 
chiral charge density is the same with the third and fourth components of the induced EM current in the leading 
contributions in 5: 
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where we have used the result from Eq. ( 26 ) . Eq. ( 27 1 is consistent with that given in Ref. [TH [5D] for a homogeneous 
magnetic field, if A3 = A4. 
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So far, we have derived the expressions for three relevant quantities using the effective action. Here is one caveat. 
In this instanton approach, the interaction of the quarks are nonlocal, resulting in that the vector and axial-vector 
currents are not conserved in a usual manner [21]. To overcome this problem, one has to gauge the quark momentum 
even in the momentum-dependent quark mass M in Eq. (Ill with the external vector and axial- vector currents. This 
treatment will provide additional terms to Eqs. (151, (19 1 and (23). We verified, however, that these terms are odd 
in the momentum integral, such as J d 4 k k^f(k 2 ), if we take only the leading contributions up to 0(5 2 ). Thus, their 
contributions can be ignored rather safely, although it is slightly problematic in the current conservation issue for all 
orders in S. 

Finally, we close this Section with remarks on the Landau levels and corresponding dimensional reduction, which 
have not been taken into account in the present work. In the presence of the strong magnetic field BqX^, the transverse 
momenta of fermions are quantized so that they are decoupled from the longitudinal components. Especially, as for 
the lowest Landau level (LLL) approximation, the transverse fermion propagator, which is factorized from that of the 
longitudinal one, can be written in a Gaussian function, 5^ LL = exp(— kj_/Bo). This makes the calculation of matrix 
elements greatly simple in such way that the integral over the transverse momenta becomes just Bo/n, where the n 
indicates the number of the fermion propagators in a matrix element |2QII22j . Hence, one is left with only the integrals 
over the longitudinal momenta. This is called the dimensional reduction in the strong magnetic field. However, our 
approach given in this Section did not show the dimensional reduction, since we limited ourselves only up to 0{5 ) in 
expanding the quark propagator and did not consider the transverse momenta quantization. The present approach 
may correspond to the effective action method given in Refs. [T21 [TBI IT9"] . 



IV. CHIRAL MAGNETIC EFFECT AT LOW TEMPERATURE 



To investigate the physical quantities in hand at low but finite temperature (T < T*), we want to discuss briefly 
how to modify the instanton variables, p and R at finite T. We will follow our previous work [T7j and Refs. [SJ [23] 
to this end. Usually, there are two different instanton configurations at finite T, being periodic in Euclidean time, 
with trivial and nontrivial holnomies. They are called the Harrington-Shepard (HS) [23] and Kraan-Baal (KB) 
calorons [23], respectively. In fact, the nontrivial holonomy can be identified as the Polyakov line as an order parameter 
for deconfincment phase of QCD. However, since we are not interested in the deconfinement phase in this work, we 
choose the HS caloron for the modification at finite T. Note that here are several caveats; 1) these modifications 
are done for a pure glue system without dynamical quarks. Hence, the instanton variables may change, if one takes 
into account the dynamical-quark contributions in the instanton distribution function. 2) Moreover, we assume the 
CP-invariant vacuum for the modifications, i.e. A = 0, whereas we are interested in the physical quantities for A^O 
for the CME. Correcting these inconsistencies may give rise to changes in the final results, although they seem small 
considering that the order of the CP- violation in reality is very tiny, but it must be beyond our scope in the present 
work, and we want to leave them for future studies. Keeping this issue in mind, we write the instanton distribution 
function at finite T with the HS caloron as follows: 

d(p, T) = C Nc A b RS N ° p b ~ 5 exp [-{A N T 2 + 7 np 2 )p 2 } . (28) 
c 

Here, the abbreviated notations are also given as: 

4 fin P -l-68aasJVc 
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Note that we defined the one-loop inverse charge and at a certain phenomenological cutoff value p cut and (R) w R. 
As will be shown, only is relevant in the following discussions and will be fixed self-consistently within the present 
framework. The Ars stands for a scale, depending on a renormalization scheme, whereas the V3 for the three- 
dimensional volume. Using the instanton distribution function in Eq. (28), we can compute the average value of the 
instanton size, p 2 straightforwardly as follows |15j : 



p 2 {T) 
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where v — {b — 4)/2. Substituting Eq. (31 I into Eq. (28 1, the distribution function can be evaluated further as: 

i 

~1 



1 



d(p, T) = C p b - 5 exp [-M(T)p 2 ] , M(T) = ^A Nc T~ 



- a2 T 4 



(32) 



The instanton number density n can be computed self-consistently as a function of T, using the following equation: 

n»M{T) = [CT{v)\» , (33) 

where we have replaced NT/V3 —> n, an d r(^ ) in dicates the T-fucntion with an argument v. Note that the C and /3 
can be determined easily using Eqs. (311 and (33), incorporating the vacuum values of the n and p: C k, 9.81 x 1CP 4 
and P ~ 9.19. At the same time, using these results, we can obtain the average instanton size p as a function of T 
with Eq. pi} . 

Finally, in order to estimate the T-dependence of the constituent-quark mass Mo, it is necessary to consider the 
normalized distribution function, defined as follows: 



d N (p,T) 



d(p,T) 
Jdpd(p,T) 



P 



h-5 



M»(T)exp [-M(T)p 2 ] 



I» 



(34) 



Now, we want to employ the large- N c limit to simplify the expression of dn{p, T). Since the parameter b is in the 
order of 0(N C ) as shown in Eq. (29 1, it becomes infinity as N c — ► 00, and the same for v. In this limit, as understood 
from Eq. ( |34[ ), d^ip^T) can be approximated as a (^-function |14j : 

lira d N (p,T) = S[ P -p(T)}. (35) 

Considering the constituent-quark mass is represented by [T3] 

Mo oc Vn / rfpp 2 5[p - p(T)] = V^(T) p 2 (T), 

we can modify the Mo as a function of T as follows: 

'VW)P 2 {T)~ 



(36) 



M n 
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7^0)^(0) 
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where we will use M « 350 MeV as done for zero T. The numerical results for the normalized p/po and n/n-o as 
functions of T in the panel of Fig. [l] As shown there, these quantities are decreasing with respect to T as expected. 
However, even beyond the critical T for the chiral restoration T c x w Aqcd ~ 200 MeV, the instanton contribution 
remains finite. In the right panel of figure, we draw the quark mass as a function of T and absolute value of three 
momentum of a quark | k \ : 



M(\k\,T) = M (T) 



2 + p 2 (T)|fc| 2 



(38) 



Note that we have ignored the Euclidean-time component of the four momentum in the M^x and Mk,T by setting 
&4 = 0. This tricky treatment simplifies the calculations in hand to a large extend, and we also verified that only 
small deviation appears in comparison to full calculations. Moreover, the p in Eq. ( |38[ ) is now a function of T as 
demonstrated by Eqs. (31) and (37l previously. As shown in the figure, the M(|fc|, T) is a smoothly decreasing 
function of T and |fc|, indicating that the effect of the instanton is diminished. For more details, one can refer our 
previous work [T?] . 

In order to compute the relevant quantities for the CME, given in the previous Section, at finite T, it is necessary to 
compute the T a and Tb in Eqs. (17) and (20 1 as functions of T. For this purpose, we employ the fermionic Matsubara 
formula with the following convention: 



/ 



d 4 fc 



/(*)-» r 



71— — OO 



d 3 k 
(2^)3 



f(w n ,k), 



(39) 



where the antiperiodic Matsubara frequency is assigned as w n = (2n + 1)ttT, where n € 1. Here, we redefine the 
mass-related functions in Eq. (14 1 for the case at finite T as follows: 



(40) 
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FIG. 1: Left: normalized p/po and n/no as functions of T for N c — 3. Right: M as a function of T and absolute value of the 
momentum fc|. 



It is worth mentioning that, in general, the Lorentz invariance is broken down due to the periodic Euclidean time 
for the direction p = 4, i.e. the Matsubara frequency. Hence, the electric part of the field strength (i 7 ^ and F^ where 
i = 1 *~ 3) for T ^ becomes different from that for 7 1 = 0. However, we are free from worrying about such a case 
in the present work, since we have — Fa = in our setting of the external EM field as shown in Eq. (16). As 
a result, the absolute values of the (j3)f.<s and (j±)F,& m Eq. (23 1 are the same even for 0, and the (JiJfs and 
O2) f,s also remain unchanged. 

Now, we are in a position to compute the T a ,b in Eqs. (17 1 and (20) as functions of T. It is convenient to employ 
a summation identity for this purpose: 



E 



AT 



1 



E 2 ) 2 



E 2 



1 (1 



-E k , T /T\ 



1 



r/T 



E k .T (1 + e-^.WT) 2T (1 



M k - 



(41) 



where we employed a notation Mk,T for the summation, and the quark energy is written as E\ T — k 2 + (1 + <5 2 )M| T . 
Then, one arrives finally at 



Fa = N c 



d 3 k 

(2tt)» 



M fe ,T I I + Mu T Mu T ) Mk, T , 3=i 



d 3 k 



-Mk.TM k TM k T- 



(42) 



By substituting these T a ,b into Eqs. (151, (19 1, and (22), we can obtain the v.e.v. of the chiral density, chiral charge 
density, and the induced EM current as functions of T. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we present the numerical results for the v.e.v. of the (p x ), (n x ), and (j^) with the nontrivial 
topological charge Q t ^ in the external magnetic field B = B x 3 . First, we present the numerical results for 

) as a function of Bq in the left panel of Fig. [2] There, we depict them 
for different temperatures, T = (0, 50, 100, 150, 200) MeV, separately. As a trial, we choose e = 10~ 3 , which gives 
8 w 0.0136 in Eq. pj. In other words, the P- and CP- violation effects are about 1% in comparison to non- violating 
quantities. Moreover, since the (p x ) is almost linearly proportional to e 2 , one can easily estimate it for different e (or 
5) values. Obviously as shown in the figure, the \(p x )\ grows rapidly with respect to Bo, manifesting the CME. As T 
increases, the strength of the curve decreases. At the same time, the slope of the curve, d\(p x )F,s\/dB , gets smaller. 
This tendency means that the \(p x )f,s\ becomes insensitive to the external magnetic field as T gets higher, and the 
same for the CME. The reason for this tendency can be understood by the decreasing instanton effect. 

In the right panel of Fig. [2j we show the \(p x )f,s\ as a function of T for different strengths of the magnetic fields, 
£?o = (0, 0.5, 1.0, 1.5, 2.0) GeV 2 . As shown there, the strength of a curve depends on that of the magnetic field. As 
T increases, the becomes small, indicating that the instanton effect is reduced. In other words, decreasing 

tunneling effect corresponding to Q t — > 0. Thus, unless there is another mechanism to make the Q t nontrivial, such 
as the sphaleron, the |(/9 x )f,<s| as well as the CME decrease monotonically with respect to T for a finite Bq value. 



the local chiral density \(p x )f.s\ in Eq. (15 
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Bo [GeV 2 ] T [GeV] 

FIG. 2: Chiral local density, \(p x ) F ,s\ x 10 4 , as functions of B (left) and T (right) for e = 1CT 3 . The five curves in each panel 
correspond to those for T = (0,50^100,150,200) MeV (left, from top to bottom) and B = (0, 0.5, 1.0, 1.5, 2.0) GeV 2 (right, 
from bottom to top). 



Interestingly, the \(p x )f,s\ decreases much faster for the larger Bq value, since the T-dependence of the quantity 
becomes more obvious and strengthened. 

If we compare the present results to a recent lattice QCD simulation given in Ref. |llj . in which the correlation of 
the chiral density (p x ) was taken into account, we observe similarity that the local chiral density becomes insensitive 
to the magnetic field as T increases: the decreasing slope. However, since we have considered only the contributions 
proportional to the Bq for the (p x ), we can not reproduce the nonzero values for it at Bq rj as shown in Ref. |llj . 
We note that, although we did not include in the present numerical calculations, there are additional terms, which 
are independent on the B in expanding Eq. ( 12 1 : 



( Px )s = AN C 



d i k 



M 



8S 2 k 2 M 



k 2 + (1 + S 2 )M 2 [k 2 + (1 + 6 2 )M 



212 



(43) 



Obviously, the first term in the square bracket in the r.h.s. of Eq. ( |43[ ) is just the chiral condensate with the nontrivial 
Qt. Note that the second term is far smaller than the first one (o~<C 1). Since the chiral condensate is a smoothly 
decreasing function of T |17j . it will give different and nonzero values for the (p y ) p x fo r each T at Bq = 0. 

Now, we are in a position to take into account the (n x ), (j^), and fi x in Eq. (26). However, as understood by 
Eqs. (23) and (27 1, these quantities are related to each other, beside the differences in the constant multiplying- 
factors in terms of the S and Bo. Hence, it is enough to consider only the (j'3,4) f,S to see the overall behaviors of 
the quantities. Since the (^'3,4)^,5 is functions of undetermined quantities which have been assumed to be small 
along the ^-direction, we compute the absolute value of the normalized quantity | (73,4)^5/(^3,4) | for convenience. 

In the left panel of Fig. |J we draw \{h A / {iA A ,?,)) f,s\ as functions of B a for T '= (0, 50, 100, 150, 200) MeV, 
separately. The curves are increasing linearly with respect to Bo as expected, indicating the CME. Being similar to 
the local chiral density, the | (j3 ; 4/(^4, 3)}f,s\ becomes insensitive to the magnetic field, as T increases. Hence, we can 
conclude that the signal of the CME is weakened with respect to T at a certain value of Bo- In Fig. |4j we draw the 
\(j3,4:/ {^4,3)} f,s\ as functions of B and T to sec its overall behavior. Although there can be a increasing competition 
between the instanton and sphaleron contributions in the high-T region (T > T c x ), it turns out that, in our present 
framework, the instanton contribution to the CME becomes about one-tenth at T k 280 MeV in comparison to its 
maximum and almost disappears above T « 300 MeV as show in the figure. 



VI. SUMMARY AND CONCLUSION 



In the present work, we have investigated the chiral magnetic effect (CME), which signals the nontrivial topological 
charge Q t of the QCD vacuum. The evidence of the CME can be explored by measuring the electromagnetic (EM) 
current induced by the external magnetic field, which is generated from the non-central heavy-ion collision. We were 
interested in the CME in the low-temperature region below the chiral phase transition, T < T* w 200 MeV. To this 
end, we employed the instanton vacuum configuration with the nonzero instanton-number fluctuation A, which is 
proportional to the Qt- In order to incorporate with finite T, we made use of the Harrington-Shepard caloron with 
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B [GeV 2 ] T [GeV] 

FIG. 3: Induced EM current indicating the CME, \(j 3 ,4,/ (1^3)) F>s \ X 10 3 , as functions of B (left) and T (right) for e = 10" 3 . 
The five curves in each panel correspond to those for T — (0, 50, 100, 150, 200) MeV (left, from top to bottom) and -Bo = 
(0, 0.5, 1.0, 1.5, 2.0) GeV 2 (right, from bottom to top). 




FIG. 4: Induced EM current indicating the CME, \(j3A/{iA 4 ^)} F ^\ x 10 3 , as functions of B and T for e = 10 3 . 

the trivial holonomy for the T-dependent constituent-quark mass. The fermionic Matsubara formula was taken into 
account for the anti-periodic sum over the Euclidean time direction. We calculated the vacuum expectation values 
of the local chiral density, chiral charge density, and EM current, which are induced by the external magnetic field 
with the nontrivial Q t , using the linear Schwingcr method and functional derivatives of the effective action with the 
external sources. We expanded the relevant matrix elements up to 0(5 2 ), where the 6 represents the strength of the 
P- and CP- violation and is proportional to the Qf We have observed the followings: 

• The transverse EM current induced by the external magnetic field is far smaller than that of the longitudinal 
one, j|| S> j±, and their ratio is proportional to 8, \j\\/j±\ oc 5 ~ Qt |11| . 

• We obtain a relation such that the strength of the chiral charge density is the same with that of the induced EM 
current in the £3- and i^-directions: (n x ) = | {^3,4} | as long as the Lorentz invariance remains unbroken |12[I13|. 

• The correct relation between the fourth component of the EM field and the the chiral chemical potential is 
found: i5A 4 = [i x [15]. 

• The CME, represented by the local chiral density and induced EM current, becomes insensitive to the external 
magnetic field as T increases. This can be understood by that the instanton, i.e. tunneling effect gets diminished 
with respect to T: Q t -> [TT]. 
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• The local chiral density and induced EM current as functions of T decrease faster for the larger external magnetic 
field, since the decreasing instanton effect is strengthened. 

• The induced EM current from the instanton contribution decreases smoothly with respect to T and becomes 
one-tenth at T w 280 MeV and negligible beyond T sa 300 MeV, if we ignore the sphaleron effect. 

To have more information on the CME and compare with the lattice simulation, we want to investigate the correlation 
functions of the chiral density and the induced EM current as done in Ref . [TT] . Related works are under progress 
and appear elsewhere. 



Acknowledgment 

The author thanks K. Fukushima, M. M. Musakhanov, and C. W. Kao for fruitful discussions. He is also grateful 
to K. S. Choi for technical supports on the numerical calculations. This work was supported by the NSC96-2112- 
M033-003-MY3 from the National Science Council (NSC) of Taiwan. 



[1] 


D. 


[2] 


T. 


[3] 


P. 


[4] 


M 


[5] 


D. 


[6] 


Y. 


m 


A. 


[8] 


S. 


[9] 


D. 


[10] 


S. 


[11] 


P. 


[12] 


K. 


[13] 


E. 


[14] 


D. 


[15] 


T. 


[16] 


H. 


[17] 


S. 


[18] 


J. 


[19] 


J. 


[20] 


N. 


[21] 


S. 


[22] 


D. 


[23] 


B. 


[24] 


T. 



Diakonov, Prog. Part. Nucl. Phys. 51, 173 (2003). 
Schafer and E. V. Shuryak, Phys. Rev. D 53, 6522 (1996). 
Arnold and L. D. McLerran, Phys. Rev. D 37, 1020 (1988). 
. Fukugita and T. Yanagida, Phys. Rev. D 42, 1285 (1990). 
Diakonov and A. D. Mirlin, Phys. Lett. B 203, 299 (1988). 
Maezawa et al, J. Phys. G 34, S651 (2007). 

Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 63, 034502 (2001). 
A. Voloshin, Phys. Rev. C 70, 057901 (2004). 

E. Kharzeev, L. D. McLerran and H. J. Warringa, N ucl. Phys. A 803, 227 (2008). 
A. Voloshin [STAR C ollaboration] , |arXiv:0806.0029| [nucl-exj . 

V. Buividovich et aL, [a"rXiv:0907.0494| [hepdatj.~ 

Fukushima, D. E. Kharzeev and H. J. Warringa, Phys. Rev. D 78, 074033 (2008). 

D'Hoker and J. Goldstone, Phys. Lett. B 158, 429 (1985). 

Diakonov, M. V. Polyakov and C. Weiss, Nucl. Phys. B 461, 539 (1996). 

Schafer and E. V. Shuryak, R ev. Mod. Phys. 70, 323 (1998). 

J. War ringa, |ajX iv:0906.2803 [hep-ph]. 

i. Nam, arXiv:0905.3609 [hep-ph]. 

S. Schwinger, Phys. Rev. 82, 664 (1951). 

F. Nieves and P. B. Pal, Phys. Rev. D 73, 105003 (2006) 
Sadooghi and A. Jafari Salim, Phys. Rev. D 74, 085032 (2006). 
i. Nam and H. -Ch. Kim, Phys. Rev. D 74, 076005 (2006) 

K. Hong, Y. Kim and S. J. Sin, Phys. Rev. D 54, 7879 (1996). 
J. Harrington and H. K. Shepard, Nucl. Phys. B 124, 409 (1977). 
C. Kraan and P. van Baal, Nucl. Phys. B 533, 627 (1998). 



